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Abstract 

In the present paper we consider the bouncing braneworld scenario, in which the 
bulk is given by a five-dimensional charged AdS black hole spacetime with matter 
field confined in a -D3 brane. Then, we study the stability of solutions with respect 
to homogeneous and isotropic perturbations. Specifically, the AdS black hole with 
zero ADM mass and charge, and open horizon is an attractor, while the charged 
AdS black hole with zero ADM mass and flat horizon, is a repeller. 
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1 Introduction 



Motivated by string/M theory, the AdS/CFT correspondence, and the hierarchy problem 
of particle physics, braneworld models were studied actively in recent years [l]-[4]. In 
these models, our universe is realized as a boundary of a higher dimensional spacetime. 
In particular, a well studied example is when the bulk is an AdS space. In the cosmological 
context, embedding of a four dimensional Friedmann- Robert son- Walker universe was also 
considered when the bulk is described by AdS or AdS black hole [5, 6]. In the latter case, 
the mass of the black hole was found to effectively act as an energy density on the brane 
with the same equation of state of radiation. Representing radiation as conformal matter 
and exploiting AdS/CFT correspondence, the Cardy-Verlinde formula [7] for the entropy 
was found for the universe ( see [8] , for the entropy formula in the case of dS black hole ) . 
In either of the above cases, however, the cosmological evolution on the brane is modified 
at small scales. In particular, if the bulk space is taken to be an AdS black hole with 
charge, the universe can 'bounce' [9]. That is, the brane makes a smooth transition 
from a contracting phase to an expanding phase. From a four-dimensional point of view, 
singularity theorems [10] suggest that such a bounce cannot occur as long as certain energy 
conditions apply. Hence, a key ingredient in producing the bounce is the fact that the 
bulk geometry may contribute a negative energy density to the effective stress-energy on 
the brane [11]. At first sight these bouncing braneworlds are quite remarkable, since they 
provide a context in which the evolution evades any cosmological singularities while the 
dynamics is still controlled by a simple (orthodox) effective action. In particular, it seems 
that one can perform reliable calculations without deliberation on the effects of quantum 
gravity or the details of the ultimate underlying theory. Hence, several authors [12, 13, 
14, 15] have pursued further developments for these bouncing braneworlds. However the 
authors of [16] have found that generically these cosmologies are in fact singular. In 
particular, they have shown that a bouncing brane must cross the Cauchy horizon in the 
bulk space. However, the latter surface is unstable when arbitrarily small excitations are 
introduced in the bulk spacetime. 

In this paper we describe solutions of the bouncing braneworld theory and also determine 
their stability. To do this, we use a set of convenient phase-space variables similar to those 
introduced in [18, 19]. The critical points of the system of differential equations in the 
space of these variables describe interesting non-static solutions. A method for evaluating 
the eigenvalues of the critical points of the Friedmann and Bianchi models was introduced 
by Goliath and Ellis [19] and further used in the analysis by Campos and Sopuerta [18] of 
the Randall Sundrum braneworld theory. These latter authors gave a complete description 
of stationary points in an appropriately chosen phase space of the cosmological setup 
and investigated their stability with respect to homogeneous and isotropic perturbations. 
The authors worked in the frames of the Randall Sundrum braneworld theory without 
the scalar-curvature term in the action for the brane. Cosmological solutions and their 
stability with respect to homogeneous and isotropic perturbations in the braneworld model 
with the scalar-curvature term in the action for the brane have further been studied by 
Iakubovskyi and Shtanov [20]. 
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2 D 3 Brane with Radiative Matter 



We start by considering a 4— dimensional brane in a space-time described by a 5— dimensional 
charged AdS black hole. The background metric is then given by [17] 

da 2 

dsl = -h{a)dt 2 + — + a 2 dQl 3 , (1) 

where 

, . , , uJiM 3u 2 Q 2 a 2 

in which M and Q are ADM mass and charge, respectively, L is the curvature radius of 
space, A; is a discrete parameter which indicates the horizon topology, namely k — +1, 
and -1 for a spherical, flat and hyperbolic geometry respectively, and 

-4 = -3^. (3) 

Here G 5 is the five- dimensional Newton constant, and V 3 is the dimensionless volume 
element associated with dVL 2 k 3 (a three-dimensional spacelike hypersurface of constant 
curvature) . 

Tuning the brane cosmological constant to zero, we get the following induced equation on 
the brane: 

2 k u 4 M 3u 2 Q 2 

where H — - is the Hubble parameter and dot denotes derivation with respect to the 
brane time r. If we consider the more physically relevant case in which a perfect fluid with 
equation of state of radiation is present on the brane, then the first Friedmann equation 
takes the following form 

where p is the energy density defined by p = E4/V and the four-dimensional energy 
is given by 

E, = = l -E, (6) 

where E is the thermodynamical energy of the black hole. Also in Eq.(5) po is the tension 
of the brane, while pb r is the energy density of the radiation. The Hubble equation can 
be rewritten as 

-k 8nG 4 A 4 Stt p 2 br 3uj 2 Q 2 
H = ^ + ^ P + Y + 3M| ( 2^ + ^ ) -T6^' (7) 

where 

47rpo 3 3 

A4 " ~Mf -T 2 - Ahr ~T 2 > (8) 

is the effective cosmological constant of the brane. By tuning the bulk cosmological 
constant and the brane tension At r , the effective four dimensional cosmological constant 



3 



A4 can be set to zero. This critical brane is what we would like to consider in the following. 
The radiation energy density pb r on the brane is 

Pbr = (9) 

where p r is a constant, hence we can rewrite the cosmological equation (7) as 



H 2_ ~k , M Snp r 1 3ujQ 2 Aix P 2 r 



By defining 

A=( W4 M+g|), (11) 
B = ^f' (12) 

and 

Eq.(lO) takes on the following form in terms of the above parameters 

tt2 -k A B C 

H " + ^-^ (14) 

3 Stability of the Bouncing Solutions 

In this section, we describe solutions of the braneworld theory under investigation and 
also determine their stability. To do this, we use a set of convenient phase-space variables 
similar to those introduced in [18, 19]. Now, we introduce the notation similar to those 
of [18] 

^ k = 2 Til = _7 2"' ^ A = 4172' ^ B = 6 112' ^ C = 8U2' 

a l H l a 1 a 4 H z a b H 2 a s H z 

and work in the 4-dimensional f2-space (f^, ^b, ^c)- in this space, the fl parameters 
are not independent since the Friedmann equation (4) reads 

n k + n A + n B + n c = i. (16) 

In this case the state space defined by the variables (f^, CIa, &b, &c) is 110 longer compact 
( because now Qb < ). However, we can introduce another set of variables describing 
a compact state space. Firstly, instead of using the Hubble function H we will use the 
following quantity 

D = V^+I' (17) 

and from it, let us define the following dimensionless variables 

Z = ^, (18) 
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«* = aW- (19 > 

ft. = ^ < 2 °) 

ftc = ^- PI) 

From these definitions we see that now the case H = is included. Moreover, the 
Friedmann equation takes the following form 

n k + n A + n c = 1, (22) 

which, together with the fact that — 1<Z<1 [see Equation (17)], implies that the state 
space defined by the new variables is indeed compact. 
Introducing the primed time derivative 



(23) 



~ Ddi' 

one obtains the system of first-order differential equations 

D' = -ZD[l + (q-2Q B )Z 2 }, 
Z'^-Z 2 [q-(q-2Q B )Z 2 }, 

n' k = 2n k (q-2n B )Z 3 , (24) 
Q' A = 20i[-l + (q~ 2Q B )Z 2 ]Z , 

n' c = 2n c [-3 + (q- 2n B )z 2 ]z, 

where 

l + (q- 2)Z 2 = AQ C - fifc. (25) 

The evolution equation for D is not coupled to the rest, so we will not consider it for 
the dynamical study. Therefore, we just study the dynamical system for the variables 
Cl = (Z, Cl k , Qa, &c), determined by the equations (24). 

The behavior of this system of equations in the neighborhood of its stationary point 
is determined by the corresponding matrix of its linearization. The real parts of its 
eigenvalues tell us whether the corresponding cosmological solution is stable or unstable 
with respect to the homogeneous perturbations [20]. To begin with, we have to find the 
critical points of this dynamical system, which can be written in vector form as follows 

Q 1 = f(Q), (26) 

where / can be extracted from(24). The critical points, Q*, namely the points at which 
the system will stay provided it is initially at there, are given by the condition 

m*) = 0. (27) 
Their dynamical character is determined by the eigenvalues of the matrix 

df 



\n=n*- 



(28) 
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If the real part of the eigenvalues of a critical point is not zero, the point is said to be 
hyperbolic [18]. In this case, the dynamical character of the critical point is determined 
by the sign of the real part of the eigenvalues: If all of them are positive, the point is said 
to be a repeller, because arbitrarily small deviations from this point will move the system 
away from this state. If all of them are negative the point is called an attractor because 
if we move the system slightly from this point in an arbitrary way, it will return to it. 
Otherwise, we say the critical point is a saddle point. 
We construct our models as follows: 

(1) The model k e , or (Z, fife, fi^, fic) = (e, 1, 0, 0), where e = sgn(Z). We have 

? = 2--^ = 0, (29) 

and the eigenvalues are 

(A z , A fc , \ A , A c ) = (0, 0, -AZ, -12Z) = -e(0, 0, 4, 12) (30) 

(2) The model A e , or (Z, fi fe , fi A , fi c ) = (e, 0, 1, 0). We have 

9 = 2-^ = 1. ( 31 ) 

and the eigenvalues are 

(X z , X k , X A , A c ) = (2Z, 2Z, 0, -AZ) = e(2, 2, 0, -4) (32) 

(3) The model C t , or (Z, fi fc , tt A , fi c ) = (e, 0, 0, 1). We have 

g = 2 + A = 5, (33) 

and the eigenvalues are 

(Xz, X k , X A , Ac) = (10Z, 10Z, 8Z, AZ) = e(10, 10, 8, 4) (34) 

(4) The model O, or 

(Z, £l k , Q A , Qc) = (0, fi^, fi^, fi c ) 

where fife,fi^ and fi c are constants satisfying (22) and (27). The eigenvalues of these 
points can be obtained in a straightforward manner and show a saddle point, which has 
not been included here. Hence, O represents a set of infinite saddle points whose line 
element is that of an open universe (k — — 1) with H = 0. 

The dynamical system (24) has three hyperbolic critical points as follows: 

i) The model k (k = -1), 

A = B = C = 0, a(t) = t, 
with the critical point of an attractor type. 

ii) The model A, 

k = B = C = 0, a(t) = (A) 1/4 v / 2t, 
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with the critical point of a saddle point type. 

iii) The model C, 

k = A = B = 0, a(t) = (AVCt) 1/4 , 

with the critical point of a repller type. 

iv) The model O, 

Z = 0, fl^, fl* A , Qq, 

with the critical point of a saddle point type. 

The four models have been depicted within the compact state space, in Fig.(l). There 
are just trajectories on the planes, which are invariant submanifolds of the state space. 

4 Conclusion 

In this work, we have studied bouncing cosmological solutions and their stability with 
respect to homogeneous and isotropic perturbations in a braneworld theory in which the 
bulk is given by a five-dimensional charged AdS black hole spacetime with matter field 
confined in a D3 brane. The effects of the charge of five-dimensional black hole in the 
bulk has been considered. By including this quantity in the analysis we have obtained 
four models with the critical points of an attractor, a couple saddle point and a repeller, 
respectively, and constructed the complete state space for these cosmological models. 
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Figure Captions 



Figure 1. State space for the bouncing braneworld models. The points k + ,A + and C + 
describe the critical points of an attractor, a saddle point and a repeller, respectively. 
The points k_,A_ and C_ describe the critical points of a repeller, a saddle point and 
an attractor, respectively. The points O represent a set of infinite saddle points. Only 
trajectories on the invariant planes, (Ct k = 0, Qa = 0, and Cl c = 0) which outline the 
whole dynamics, are drawn. 
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This figure "bouncing.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/hep-th/0605081v3 



